Shape of phospholipid/surfactant mixed micelles: cylinders or discs? 

Theoretical analysis. 



M.M.Kozlov 1 , D.Lichtenberg 1 and D.Andelman 2 



Department of Physiology and Pharmacology, Sackler School of Medicine, Tel Aviv 
University, Ramat-Aviv 69978, Israel 

2 School of Physics and Astronomy, Raymond and Beverly Sackler Faculty of Exact Sci- 
ences, Tel Aviv University, Ramat-Aviv 69978, Israel 



Abstract. 

We develop a theoretical model for the solubilization of phospholipid bilayers by micelle- 
forming surfactants. Cylindrical micelles, disc-like micelles and spherical micelles are 
considered as alternative resultant structures. The main question addressed is what kind 
of micelles can be expected under various thermodynamical conditions. Our analysis is 
based on a theoretical model that accounts for Helfrich energy of curvature of amphiphile 
monolayers and for the entropy of mixing of lipids and surfactants in mixed aggregates. 

We conclude that for usual values of the elastic parameters of amphiphile monolayers 
cylindrical micelles are the most probable aggregates resulting from micellization of phos- 
pholipid by surfactants. This conclusion is consistent with available experimental data. 
Conditions of formation of disc-like and spherical micelles are also determined. 
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Introduction. 



Amphiphiles tend to self assemble in aqueous solutions, mostly due to the hydrophobic 
effect 1 . In the resultant aggregates, the amphiphilic molecules are packed as monolayers, 
where their hydrophobic moieties are shielded from contact with the external aqueous 
medium by the polar head groups. Depending on their molecular structure and inter- 
actions, the amphiphiles form aggregates of different shapes . Most of the biological 
amphiphiles (phospholipids) self assemble in nearly flat bilayer membranes, forming closed 
vesicles (liposomes). By contrast, most of the commonly used surfactants form micelles 
whose radii of curvature are close to length of hydrocarbon chains 5 . 

While each of the pure compounds form in dilute solution aggregates of a particular 
type, mixtures of lipids and surfactants self assemble in either mixed liposomes or mixed 
micelles, depending on the composition 5 '. Transition from mixed bilayers to mixed 
micelles upon addition of surfactant to phospholipid vesicles, is commonly denoted as 
solubilization of the liposomes. The resultant micelles were previously described as having 
either disc-like (oblate ellipsoidal 8 ) or cylindrical shapes, 5 in apparent agreement with 
dynamic light scattering data. 

More recently, however, cryo-transmition electron microscopy, ^' 1( ^ size-exclusion high 
performance liquid chromatography, 11 small angle neutron scattering, 12 and re-evaluation 
of dynamic light scattering data 13 indicated that in most cases solubilization of liposomes 
results in formation of thread-like rather than disc-like micelles. 

The model of disc-like micelles was supported by the idea that the surfactant molecules 
form the rims of the micelles where the amphiphile monolayers are strongly curved, while 
the lipid molecules remain in flat parts of the discs. However, this qualitative consideration 
did not account for the entropy of mixing of the two components in the micelle, which 
tends to distribute uniformly the molecules of the two components over the whole surface 
of each micelle. The result of competition of these two tendencies is not obvious and 
requires a detailed theoretical analysis. 

Based on the more recent experimental results, theoretical approaches that were de- 
veloped to describe the energetics and size distributions of mixed amphiphilic aggregates 
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' b and to interpret the phase diagrams of lipid-surfactant mixtures, assumed that 
cylindrical and spherical micelles are the only possible aggregates resulting from solubi- 
lization. The question remains open whether disc-like micelles are indeed less favorable 
energetically than the cylindrical and spherical ones and if there are conditions where 
solubilization can still result in formation of disc-like micelles. 

The present work analyses and compares the conditions of surfactant-induced phase 
transition of lipid bilayers into discoidal, cylindrical and spherical micelles. Our model, 
unlike alternative theoretical approaches 16 , do not consider the detailed distribution of 
microscopic interactions in amphiphile monolayers. Instead, we describe a monolayer 
of amphiphiles by the more macroscopic Helfrich elastic model and account for the 
lipid/surfactant entropy of mixing. We show that the shape of the micelles formed upon 
solubilization of liposomes is determined by a unique parameter that depends on the tem- 
perature and on the elastic characteristics of the monolayer as expressed by its bending 
rigidity, the spontaneous curvatures of the two compounds and the Gaussian curvature 
modulus. 

Based on our calculations, we conclude that at all reasonable values of this parameter 
the predicted shapes of the mixed micelles are those of long cylinders, in agreement with 
the recent experimental data^' 1 ^ 1 . Disc-like micelles can only be expected for compounds 
whose Gaussian curvature modulus is of unusually high negative value. Another possi- 
bility to obtain disc-like micelles is to suppress the effects of the entropy of mixing by 
decreasing the temperature. 



The model. 

We consider a ternary system of water, lipid and surfactant. The concentrations of lipid 
and surfactant in water, denoted by and Np, respectively, are assumed to be much 
higher than the critical micelle concentrations (cmc). 

We consider the following states of the aqueous solution of the amphiphiles in terms 
of the most common shapes of the aggregates: flat bilayers (liposomes, whose radius is 
taken as very large in comparison to the bilayers thickness), cylindrical micelles, disc-like 
micelles, spherical micelles and mixtures of coexisting aggregates of these types. Since the 
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aim of our work is to analyze qualitatively the main pathways of solubilization of bilayers 
we do not consider the more complicated architecture of intermediate aggregates 17 such 
as ellipsoidal micelles, hyperbolic and mesh structures 17-20 . 

To characterize the composition of the system we use the area fraction occupied by 
surfactant 



where ajy and a £ are the molecular areas of the surfactant and lipid, respectively, at the 
monolayer plane. It is assumed that those specific areas do not differ for aggregates of 
various types. 

In our model, the pure lipid system (0 = 0) preferentially forms liposomes, while 
addition of the surfactant results in transition of liposomes into micelles of one of the 
shapes mentioned above. To analyze the resulting structures, we determine for each of 
the possible aggregates mentioned above the free energy as a function of composition. 
Comparing those free energies, we find for each composition the state of the lowest free 
energy, i.e. the equilibrium structure. 

The obvious difference between the structures of aggregates is the curvature of their 
monolayers. Therefore, as a first contribution to the free energy (per unit area) of the 

9 1 

monolayer, u, we consider the Helfrich energy of bending u^. The other major contribu- 
tion is that of the entropy of mixing of the two components in the monolayer, s.^ Hence, 



We neglect the translational entropy of aggregates and the entropy of polydispersity of the 
micellar sizes 14 since the related contribution to the free energy can be shown to correct 
only slightly the criteria of solubilization of the bilayers. 

The energy of bending (per unit area) of the monolayer is 



= 



a D N D 



(1) 



a D N D + a L N L 



u — ujj — T-s 




(2) 
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where c m and c p are the principal curvatures of the bent surface . The material prop- 
erties of the monolayer that determine the bending energy (2) are the bending rigidity k, 
the Gaussian curvature modulus R and the spontaneous curvature cq. 

The model (2) has been originally formulated for membrane shapes that deviates 
only slightly from flat surface . 21 Extension of this model to the cases of strongly curved 
micelles can be justified by recent analysis of the elastic properties of monolayers of 
inverted hexagonal phases, whose curvature is comparable (but opposite in sign) to the 
curvature of cylindrical micelles. 22 ' 23 

For the sake of simplicity we will assume that the moduli k and R do not depend 
on composition. While several theoretical models predict different characters of such 
dependencies, 24 ' 25 a more recent numerical calculation 26 revealed a weak dependence of 
k on the composition of a mixed membrane. Moreover, an experimental determination 
of the bending rigidity of mixed monolayers of Hjj phases did not show any pronounced 
changes of k with composition. 23 

By contrast, the spontaneous curvature cq depends strongly on the monolayer com- 
position. Although this dependence may be rather complicated, 15 ' 27 we will assume that 
co is the area-weighted average of the spontaneous curvatures of pure lipid and pure 
surfactant cjj, 

co = (1 - 4>)cl + <Pcd (3) 

This assumption is supported by numerical calculations 26 that showed that the sponta- 
neous curvature of a monolayer consisting of surfactants with different chain lengths is a 
linear function of <j) over a wide range of compositions. Furthermore, for mixed monolayers 
of H/j phases it has been shown experimentally 23 that cq depends linearly on <p. 

To account for the entropy of mixing (per unit area) of the two components we will 
use an approximate expression for ideal mixing 6 

J_ = - — [01og0 + (1 - 0) log(l - 0)] (4) 
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The following analysis depends mainly on the behavior of the entropy s near its minimal 
value at (f) — \. Within this range (4) can be substituted by a simpler equation 

J-. s = -±[20(0-1) -0.2] (5) 
k B a L 

Eq. (5) in fact approximates (4) with a good accuracy for any value of within the range 
0.2 < < 0.8, which covers the whole range of " solubilizing ratios" in any lipid/surfactant 
mixture studied thus far. 

We will assume below that the characteristic radius of the curved parts of amphiphile 
monolayers is equal p for all kinds of aggregates. It is convenient to express all the variables 
in dimensionless form. The dimensionless free energy per unit area /, the dimensionless 
curvature C, and the dimensionless temperature rj are defined as 

/ = £■«. C = pc,and, = ^. 

The dimensionless energy per unit area of the monolayer is 

/ = (Cm + Cp - Co) 2 + 2 • - • Cm • C P + V ■ <t>{4> - 1) + const (6) 

K 

where Co depends on according to (3). 

Equations (l)-(6) are the basis for the determination of the free energies of different 
micellar states and liposomes (flat bilayers) as functions of their compositions. 

Free energies. 

Disc-like micelles. 

The shape of a disc-like micelle is assumed to consist of a flat central part with radius R 
and a rim formed by a strongly curved monolayer, whose meridional principal curvature is 
1/p (Fig.l). The areas of the rim and the flat part are denoted as A e and Aj, respectively. 
The geometrical characteristics of the monolayer forming the rim are considered in details 
in Appendix A. 
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To calculate the free energy of a disc-like micelle fp as a function of its composition 
we take into account the partitioning of the surfactant between the flat part and the 
rim (Appendix B) while determining the area- weighted average of the energy (6). Denot- 
ing the dimensionless total curvature in the rim as J = (c m + c p )p, we obtain for the free 
energy (per unit area) 



f D = [(AC) 2 + v]<t>l + (2ClAC - v - 2^AC < J >)<t>D 



At 



<J 2 >+ (A0 .,o <J> 2 



V + (AC) 2 V + (AC) 2 



A? 77 + (AC) 2 A t SL SL K A t 



where < J > and < J 2 > are the values of the total curvature and its square, averaged 
over the curved area forming the rim, and AC = Cd ~ (l is the difference between the 
spontaneous curvatures of surfactant and lipid expressed in dimensionless form. 



Cylindrical micelles. 

We assume that all the cylindrical micelles have the same shape consisting of a cylindrical 
part of a length I with two semi-spherical caps at the ends (Fig. 2). The radius of the cross 
section of the cylinder and of the hemispheres is assumed to be equal to p. The areas of 
the cylindrical part and the caps will be denoted as A\ and A^, respectively, and the total 
area is At = A[ + A^. 

To derive the energy of the phase of cylindrical micelles (Appendix C) we first de- 
termine the partitioning of the surfactant between the less curved cylindrical body and 
more curved semi-spherical caps. Subsequently, we average the energy (6) over the area 
of one cylindrical micelle, including its body and two caps. As a result we obtain the free 
energy per unit area as a function of the micellar composition C 

f e = [(AC) 2 + V }4 + (2ClAC - V ~ 2AC)0 C - 2AC0 C ^ 
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Spherical micelles. 

The energy per unit area of spherical micelles of radius p derived from (6) as a function 
of their composition <fi s is 



f s = [(AC) 2 + V ]<P 2 S + (2ClAC - v ~ 4AC)0 S + 4 - 4C L + 2- + Cl 



(9) 



Flat bilayer of a liposome. 

The energy per unit area of a monolayer of a flat bilayer fa as a function of its composition 
(fib is given by 



Liposomes to Micelles Transition. 

We analyze the transition from liposomes to micelles by Gibbs graphical method, illus- 
trated in Fig. 3 for one particular set of parameters of the system. The energies of the pure 
phases (7)-(10) are convex functions of composition as presented by the curves b, c and d 
for the bilayers (liposomes), cylindrical micelles and disc-like micelles, respectively. The 
energies of mixtures of coexisting phases are presented by common tangents (dashed lines 
on Fig. 3). Compositions (fi* at the points where the common tangents touch the energy 
of bilayers (curve b) indicate phase transitions. We denote the compositions determining 
the transitions of liposomes to cylindrical, disc-like and spherical micelles, by <fi*{b — > c), 
<fi*(b — > d) and <fi*(b — > s) respectively. 

The micelles for which the composition of transition (fi* has the lowest value are 
expected to be formed upon solubilization. 

The details on determination of the compositions of transition of the basis of Eqs.(7)- 
(10) are described in Appendix D. 



h = [(AC) 2 + v}4 + (2ClAC - v)h + Cl 



(10) 
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We obtain for the transitions from bilayers to cylindrical micelles: 



V + (AC) 1 



1 



| R 77 + (AQ 2 47T P 2 



2 A^ + A; K 7] Aft + ^ 



for the transitions from bilayers to disc-like micelles: 



V + (AC) 1 



a + Ac-0*(6-d) 



1 <J 2 > | g ^ + (AQ 2 4tt P 2 

2 < J > K 7] A e < J > ^ ' 



and for the transition from bilayers to spherical micelles: 



77 



Ci + AC-0*(6^s) 



1 | 1« 77 + (AQ 2 
2 k 77 



(13) 



Criteria for the shape of micelles. 



To compare the critical compositions given by Eqs.(ll)-(13), we use the expressions for 
A e , < J 2 > and < J > derived in Appendix A. 

It follows from (11)-(13) that the type of the formed micellar phase is controlled by 
a unique parameter 

A^." +(A ° 2 (14) 

K 7/ 

To illustrate it, we show in FigA (A to C) the dependence of the critical composition for 
the transitions of the bilayer into cylindrical micelles <f)*(b — > c) and into disc-like micelles 
4>*(b — > d) as functions of the micellar area At for different values of A (A = —0.3 in panel 
A, A = —0.6 in panel B, A = —2. in panel C). 

For comparison, we show on the same figure the critical composition for the transition 
into spherical micelles 0* {b — > s) presented as a constant dashed line as the size of micelles 
of this type is fixed. As obvious from Figs. 4, for A = —0.3 (Fig.4A) the lowest critical 
composition corresponds to the transition into long cylindrical micelles; for A = —0.6 
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(Fig.4B) the lowest critical composition corresponding to the transitions into the disc- 
like micelles of finite area; and, finally, in the case A = —2. (Fig.4C) the lowest critical 
composition is equal zero and corresponds to the transition into spherical micelles. 

Detailed analysis of the model shows that there are two values of the parameter A, 
equal to Ai = —1/2 and A2 = —2 that separate the solubilization of bilayer into three 
different regimes. For 

A>-l/2 (15) 



the phase transition results in formation of long cylindrical micelles; 
for 



-2 < A < -1/2 (16) 



the transition leads to formation of disc-like micelles of a finite radius, 
while for 



A < -2 (17) 



the bilayer transforms into spherical micelles. This is illustrated in Fig. 5, which shows 
the type of phase transition for different values of the parameter A. 

Discussion. 

We have shown that the type of micelles resulting from solubilization of bilayers is deter- 
mined by the value of the parameter A (14). This parameter depends on the difference of 
spontaneous curvatures of surfactant and lipid Ac = — c^, the bending rigidity k, the 
Gaussian curvature modulus R, and the temperature T 

R k B T + \k • a L • (Ac) 2 

k ' k B T 1 ' 
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To understand qualitatively these results, let us recall that according to the model of 
membrane elasticity 21 the Gaussian curvature modulus R, which controls the tendency 
of the monolayer to change its topology, can be either positive or negative. Negative 
values of R favor the division of each closed monolayer into as large as possible number of 
separated closed monolayers, whereas positive values of R result in the opposite tendency, 
i.e. in recombination of separated membranes into a single one. Hence, negative values 
of R (yielding negative values of A) prefer large number of relatively small disc-like or 
spherical micelles with respect to a fewer long cylindrical micelles. An opposite tendency 
competing with the effects of the negative R relates to the energy of bending. This 
energy is controlled by the bending modulus k, the difference in spontaneous curvatures 
of the components Ac and the effectiveness of repartitioning of the surfactant between 
the parts of the monolayers with different curvature, which, in turn is determined by the 
temperature T. The competition between these tendencies is expressed by the parameter 
A (18) and criteria (15)-(17). 

Most of the parameters needed to estimate A are known or can be estimated with 
reasonable accuracy. We will assume an area per lipid molecule = 0.6nm 2 , a radius 
of curvature of micelles p = 1.5nm, a bending modulus of the monolayer k = lOkgT 
(at room temperatures), a spontaneous curvature of lipid C£ = 0, and a spontaneous 
curvature of the surfactant cp = 1/ p. 

In the lack of a reliable experimental value for the modulus of Gaussian curvature 
of the monolayer R, we relate to the theoretical prediction of R, as derived from recently 
developed models 24 ' 28 which predict R to be negative and quite small in its absolute value. 

With these estimates we can re-express criteria (15)-(17) in terms of the value of the 
Gaussian curvature modulus R. In particular, formation of disc-like micelles occurs only 
when 

- < -0.2 (19) 

K 

Numerical calculations 24 performed on lipids and surfactants with usual characteristics 
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give larger values of R than required by (19). Therefore, formation of disc-like micelles 
seems to be a very rare event occurring only for lipid/surfactant mixtures with unusual 
properties. This prediction is in agreement with recent experimental results. However, 
we note that molecules with highly negative R, which satisfy (19), can exist. For such 
compounds, formation of disc-like micelles should be expected. 29 

In conclusion, present model provides the basis for understanding the surfactant- 
induced transformation of bilayers into various types of mixed micelles. Nonetheless, 
predictions based on this model should be considered as qualitative rather than quan- 
titative ones. Further development of the model requires more detailed experimental 
information on the elastic properties of mixed amphiphilic monolayers. 

Acknowledgment. We would like to thank R.Granek, W.Helfrich and S.Safran for 
stimulating discussions. Support from the German-Israeli Foundation (GIF) under grant 
No. 1-0197 and the Deutsche Forschungsgemeinschaft through SFB 312 is gratefully ac- 
knowledged. 



- 12 - 



References. 

1 Tanford C, The Hydrophobic Effect - Formation of Micelles and Biological Mem- 
branes 2nd.ed., John Wiley and Sons, 1980. 

2 Israelachvili, J.N.; Mitchell, D. J.; Ninham,B.W. J. Chem.Soc. Faraday Trans. II 1976, 
72, 1521. 

3 Israelachvili,J.N.; Marcelja,S.; Horn,R.G. Q. Rev. Biophys. 1980, 13, 121. 

4 Gelbart W.M.; Roux,D.; and Ben-Shaul,A. (Editors). Modern Ideas and Problems in 
Amphiphilic Sciences. Springer, 1993. 

5 Lichtenberg D. Liposomes as a model for solubilization and reconstitution of mem- 
branes. In: Handbook of Nonmedical Applications of Liposomes, v. II, Y.Barenholz 
and D.D.Lasic eds., CRC Press, Boca Raton New York London Tokyo, 1996. 

6 Andelman D.; Kozlov,M.M.; Helfrich,W. Europhys.Lett.1994, 25, 231. 

7 Fattal,D.R.; Andelman, D.; Ben-Shaul,A. Langmuir. 1995, 11, 1154. 

8 Schurtenberger,P.; Mazer,N.A.; Kanzig,W. J. Phys.Chem. 1985 89, 1042. 

9 Walter,A.; Vinson,P.K.; Kaplun,A.; Talmon,Y. Biophys. J. 1991, 60, 1315. 

10 Vinson,P.K.; TalmonY.; Walter,A. Biophys. J. 1989, 56, 669. 
Edwards, K; Almgren,M.; Bellare,J.; Brown,W. Langmuir 1989, 5, 473. 
Edwards, K; Almgren,M. J. Colloid Interface. Sci. 1991, 147, 1 

11 Nichols,J.W.; Ozarowski , J . Biochemistry 1990, 29, 4600. 

12 Pedersen,J.S.; Egelhaaf,S.U.; Schurtenberger,P. J. Phys.Chem. 1995, 99, 1299 
Hjelm,R.R; Alkan,M.H.; Thiyagarajan,P. Mol.Cryst.Liq.Cryst. 1990, 180A, 155. 

13 EgelhaaLS.U.; Schurtenberger,P. J. Phys.Chem. 1994, 98, 8560 

14 Ben-Shaul,A.; Rorman,D.H.; Hartland,G.V.; Gelbart,W.M. J. Phys.Chem. 1986, 90, 
5277. 

15 Szleifer,L; Ben-Shaul,A.; Gelbart,W.M. J.Chem.Phys. 1987, Stf, 7094. 

16 Zoeller,N.J; Blankschtein,D. Ind.Eng.Chem.Res. 1995, 34, 4150. 

17 Kekicheff,P; TiddyG.J.T. J. Phys.Chem. 1989, 55, 2520. 

18 Hyde,S.T. Pure and Appl. Chem. 1992, 64, 1617. 

19 Hyde,S.T. Colloque de Physique 1990, C7, 209. 

20 Fredrickson,G.H. Macromolecules 1991, 24, 3456. 

21 Helfrich,W. Z.Naturforsch. 1973, 2S(cj, 693. 

- 13 - 



Kozlov,M.M.; Leikin,S.; Rand,R.P. Bioph. J. 1994, 67, 1603. 
Leikin,S.; Kozlov,M.M.; Fuller,N.L.; Rand,R.P. 1996. Biophys.J. 1996, 71, 2623. 
Szleifer,!; Kramer, D.; Ben-Shaul,A.; Gelbart, W.M.; Safran,S.A. J.Chem.Phys. 
1990, 92, 6800. 

Dan,N.; Safran,S.A. Europhys.Lett. 1993, 21, 975. 
Dan,N.; Safran,S. A. Macromolecules 1994, 27, 5766. 
May,S.; Ben-Shaul,A. J.Phys.Chem. 1995, 103, 3839. 
Kozlov,M.M.; Helfrich,W. Langmmr 1992, 5, 2792. 

Safran,S.A. Statistical Thermodynamics of Surfaces, Interfaces, and Membranes. 
Addison- Wesley Publ.Comp., 1994. 

Granek,R.; Gelbart,W.M.; Bohbot,Y.; Ben-Shaul,A. J.Chem.Phys. 1994, 101, 4331. 



- 14 - 



Appendix A: Geometrical model of disc-like micelle. 

We model the rim of a disc-like micelle as a surface of revolution of a semi-circle about a 
vertical axis as illustrated on Fig.l. 

The total curvature of the rim c m + c p depends on the position along the surface de- 
termined by the angle a (Fig.l). Indeed, one principal curvature (meridional curvature) 
is constant, 

c m = - (Al) 
P 



while the second one (parallel curvature) is given by 



1 cos a . . . 

c p = 7TT ( A2 ) 

p cos a + R/p 



and changes on the rim. 

The resulting dimensionless total curvature is 



cos a + Rj p 



The element of the area of the rim is 



dA e = 2-Kp 2 ^— + cos c^j da (AA) 



and the total area of a circular rim is 



A e = Anp 2 {l + l-) {Ah) 
z p 
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Averaging the curvature and its square over the area of the rim gives 



< J >= \~f (AG) 
? + H 



<J Z >= 



9 P 

7T , P 

2 + H 



2 + 



-.tg- 



_ p - 

R 



(.47) 



Appendix B: Energy of phase of disc-like micelles. 

We consider the free energy of disc-like micelles. The shape of one micelle is illustrated on 
Fig.l. The area of monolayer forming the rim of a micelle and that forming the flat central 
part are denoted by A e and Aj, respectively. The composition of a micelle averaged over 
its entire area is 

The free energy of a micelle is equal to the sum of the free energies of the rim and of 
the flat part. Let us consider them separately. 

The free energy per unit area of the rim given by (6) and accounting for (3) and (4) 
can be written in dimensionless form as 

7=(</-Ci-0-AC) 2 + #(0-i) (si) 



where J and <fi are, respectively, the local values of the curvature and composition. We 
will add later the Gaussian curvature term as its integral over the micellar closed surface 
is equal to 8tck/ k. 

The composition of the rim averaged over its area will be denoted as < (f) >= <p e . 
Minimizing the energy (Al) at fixed average composition e we find the distribution of 



- 16 - 



the composition along the surface of the rim 



Accordingly, 



<0 2 >=0 e 2 + pPL p (<J 2 >-<J> 2 ) (53) 



and 



< J(p >= 4> e < J > + ^P +r) [< J 2 >~<J> 2 ) (BA) 



Averaging the free energy (Bl) over the area of the rim and accounting for the equations 
(B2)-(B4) we obtain 

fe = C 2 L + <J 2 >+ [(AC) 2 + ?# 2 - 2(C L + AC-0 e )< J> - 
(AC) 2 



(AC) 2 + V 



<J Z >-<J> Z + (2C X AC - rj)<t> e (B5) 



The free energy per unit area for the flat central part of a micelle (6) in a dimensionless 
form is 

f f = [(AC) 2 + v]<f>} + (2C L AC - V )<f>f + Cl (B6) 



where 0/ is the composition of the flat part. 
The total energy of a micelle is 

F D = f e A e + f f A f (B7) 
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while the compositions of the rim and the flat part satisfy the condition 

4> e A e + 4>fA f = 4> D A t (B8) 

where At = A e + Aj is the total area of monolayer forming the micelle. Minimizing the 
free energy (B7) and taking into account (B8) we obtain 

and 

The resulting expression for the free energy of a micelle is 

F d = A t -f D (511) 



where fjj is the free energy of disc-like micelles per unit area given by (7). 

The expressions for < J >, < J 2 >, A e and At for a specific model of micelle shape 
used here are derived in the Appendix A. 

Appendix C: Energy of phase of cylindrical micelles. 

Shape of a cylindrical micelle is assumed to consist of cylindrical body of length I and 
area A\ and two spherical caps of total area A^. The radius of the cylinder and the caps 
are assumed to be equal to p. Here, their area is equal to A[ = 2npl and A^ = 4-7rp 2 , 
respectively (Fig. 2). 

Composition of a cylindrical micelle averaged over its area will be denoted as (f> c , 
while 4>h an d 4>l are the compositions of the caps and cylindrical body of a micelle. In 
analogy with the calculation performed in Appendix B, the dimensionless free energy of 
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a micelle is 



fc = A h 



(2-Ci-AC-0/J 2 +#!-r^ 



+Ai 



+8tt-p 2 (CI) 

K 



The compositions e and 0; are related to the average composition <p c by 



<l>h A h + <t>l A l = <t>cM 



(C2) 



Minimization of the energy (CI) while accounting for (C2) results in expressions for the 
compositions of the caps and the cylindrical part 



( Ph = ( Pc+ / A >X 2 I T 



(C3) 



0/ = 0C ~ 



AC Aft 
(AC) 2 +77 A t 



(C4) 



Inserting (C3) and (C4) into (CI) we get the dimensionless free energy of a cylinder mi- 
celle 



F c = A t f c 



(C5) 



where / c given by (8) is the free energy per unit area of surface of cylindrical micelles. 

In analogy to the derivations above we obtain the expressions (9) and (10) for the 
free energies of spherical micelles and flat bilayers. 



Appendix D: Derivation of compositions of phase transition 
by common tangent construction. 
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The compositions of bilayers determining their transition to micelles satisfy the condi 
tions of equal and equal / — 
to determine these compositions are 



d f df 

tions of equal ^ and equal / — 0^ in the both phases. Hence, the expression necessary 



|^ = 2[(AC) 2 + r,\<f> D + (2ClAC - V ~ 2^AC < J >) 



(Dl) 



df D 
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Figure captions. 



Fig.l. Schematic representation a disc- like micelle. c m = ^ and c p are the meridional 
and parallel curvatures of the surface of the rim, respectively, and R is the radius of the 
flat central part of the micelle. Meridional angle a determines position along the profile 
of the rim. 

Fig. 2. Schematic representation of a cylindrical micelle. / is the length of cylindrical 
body, p is the radius of curvature of the semi-spherical caps. 

Fig. 3. Free energies of the different states of the system. 

- liposomes (bilayers) (curve b, according to Eq.10), 

- disc-like micelles (curve d, according to Eq.7), 

- cylindrical micelles (curve c, according to Eq.8), 

- the phase of coexisting liposomes and disc-like micelles (common tangent b-d), 

- the phase of coexisting liposomes and cylindrical micelles (common tangent b-c) 

The chosen parameters on the curves are: r] = 0.75, A( = 1., (i = 0, R = 0, the radius 
of a disc-like micelle R = 100 • p, the length of a cylindrical micelle / = 10000 • p. 
Points (1) and (2) on the 0-axis indicate the critical compositions 0*(6 — > d) and 
0*(6 — > c), respectively. 

Fig. 4. Critical compositions of transition of liposomes into disc-like micelles ( (f)*(b — ► d), 
curve d according to Eq.12) and into cylindrical micelles ( (f)*(b — > c), curve c according 
to Eq.ll) as functions of surface area of one micelle for different values of the parameter 
A. (For comparison, the dashed line s (Eq.13) shows the constant critical composition of 
transition of liposomes to spherical micelles (4>*(b — ► s)). 
(A) A = -0.3; (B) A = -0.6; (C) A = -2.; 

Fig. 5 Phase diagram of the shapes of micelles resulting from solubilization of bilayers 
(liposomes) at different values of parameter A. 
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